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A simple local proof of Noether's Second Theorem is given. This proof immedi- 
ately leads to a generalization of the theorem, yielding conservation laws and/or 
, explicit relationships between the Euler-Lagrange equations of any variational prob- 

lem whose symmetries depend upon a set of free or partly-constrained functions. 
Our approach extends further to deal with finite difference systems. The results 
are easy to apply; several well-known continuous and discrete systems are used as 
illustrations. 
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1. Introduction 



Variational methods are endemic in physics and engineering. Commonly, physical 
symmetries are built into the Lagrangian; these correspond to conservation laws via 
Noether's Theorem. In classical mechanics, for instance, invariance under transla- 
tion in time is associated with conservation of energy, invariance under rotations in 
the base space yields conservation of angular momentum, and so on. Noether's The- 
orem applies only to systems of Euler-Lagrange equations that are in Kovalevskaya 
form (see Olver (1993) for details). For other Euler-Lagrange systems, each nontriv- 
ial variational symmetry still yields a conservation law, but there is no guarantee 
£f~j ! that it will be nontrivial. 

■ It is still not widely known that Noether's famous paper on variational problems 

(Noether, 1918) includes a second theorem, which applies to variational symmetries 
that depend on arbitrary locally smooth functions of all of the independent vari- 
ables. (Henceforth, we restrict attention to functions that are locally smooth in each 
continuous variable.) Noether's Second Theorem states that such symmetries exist if 
and only if there exist differential relations between the Euler-Lagrange equations. 
In this case the Lie algebra of variational symmetries is infinite-dimensional and 
the Euler-Lagrange equations are not in Kovalevskaya form. Brading (2002) is an 
interesting and accessible history of the mathematics and physics of Noether's two 
theorems and Weyl's gauge symmetries. Brading observes that Noether's Second 
Theorem has received scant attention in the physics literature, despite its impor- 
tance for equations that have gauge symmetries. 

Noether's Theorem deals primarily with finite-dimensional Lie algebras of vari- 
ational symmetry generators, whereas Noether's Second Theorem addresses the 
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infinite-dimensional case when generators depend on arbitrary functions of all in- 
dependent variables. However, there is an intermediate case that occurs commonly: 
generators may depend upon functions that satisfy some constraints. For instance, 
the variational problem 



for the one-dimensional wave equation is invariant under Lie point transforma- 
tions whose characteristic is of the form g(x — t), where g is an arbitrary function. 
Therefore the Lie algebra of variational symmetry generators is infinite-dimensional. 
However, the differential constraint 



applies; there are no generators that depend upon an unconstrained arbitrary func- 
tion g(x,t). Rosenhaus (2002) examines this type of problem from the global view- 
point, determining which spatial boundary conditions are necessary for the existence 
of integral quantities that are conserved in time. However, the effect of constraints 
on local conservation laws has remained an open problem until now. 

Some important variational problems have variational symmetry generators that 
depend on several arbitrary functions which are linked by constraints. Noethcr's 
Second Theorem has not yet been generalized to deal with this type of problem. 

A discrete variational problem may be a model of a discrete process, or it may 
be a discretization of a continuous variational problem. Either way, discrete ana- 
logues of Noether's two theorems are important. Noether's Theorem for finite differ- 
ence variational problems has now been well-studied; see, for example, Dorodnitsyn 
(2001), Hickman & Hereman (2003) and Hydon & Mansfield (2004). If a discretiza- 
tion preserves a variational symmetry, the discrete version of Noether's Theorem 
guarantees that the corresponding conservation law is also preserved. 

Christiansen & Halvorsen (2011) have constructed a gauge symmetry-preserving 
discretization of the Maxwell-Klein-Gordon equations for a charged particle in an 
electromagnetic field. They used the discrete analogue of Noether's Theorem to 
obtain a linear relationship between the discrete Euler-Lagrange equations; this 
was viewed (correctly) as the satisfaction of a physically-important constraint. So 
there is evidently a need for a discrete analogue of Noether's Second Theorem. 

The current paper has three purposes. First, we give a simple local proof of 
Noether's Second Theorem for a single arbitrary function and extend this to deal 
with variational symmetries that depend on several arbitrary functions. Then we 
show what happens locally when the arbitrary functions are constrained in some 
way. Finally, we transfer these ideas to discrete variational problems. Some well- 
known systems are used to illustrate the application of our results. 

2. Symmetries of Lagrangians: the continuous case 

We begin by summarizing a few key facts about variational symmetries; for details, 
see Olver (1993). We work locally, using independent variables x = (x 1 ,...,x p ) 
and dependent variables u = (it 1 , ... , u q ). Derivatives of each u a arc written in the 
form u°j, where J = (J 1 ,... ,j p ) is a multi- index; here each f is a non- negative 
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integer that denotes the number of derivatives with respect to x l , and so u a = u a . 
The total derivative with respect to x l is the operator 

dx l ' dx % du", ' 
and we shall use the following shorthand: 

Dj=D(D( ...D£, (Dj)t = (-D)j = (-I)'' if Dj, 

Here f denotes the formal adjoint. Henceforth, we adopt the Einstein summation 
convention to avoid a plethora of summation symbols. Consider the action 



£[u] = J L(x, [u])dx, 



where [u] in the integrand is shorthand for u and finitely many of its derivatives. 
The Euler-Lagrange equations are obtained by varying each u a and integrating by 
parts: 

= S u jC[u] = J ^ft#dx= J -^T>j(Su a )dx = jE Q (L)6u a dx. (2.1) 

Here the Euler operator E a , which corresponds to variations in u a , is defined by 

' dL \ 



E Q (L) = (-D)j 



du a 3 



Following Olver, we have assumed that variations are restricted to ensure that there 
are no contributions from the boundary of the domain. With this formal approach, 
two Lagrangians L and L' lead to the same set of Euler-Lagrange equations, 

E o (L)=0, a=l,...,q, (2.2) 

if and only if L — V is a total divergence. 

Typically, the set of all generalized (or Lie-Backlund) symmetries of the Euler- 
Lagrange equations (|2.2[) is not a group. However, just as for point symmetries, one 
can find infinitesimal generators (in evolutionary form, prolonged to all orders), 

d d d 

x = Q a in: + (AQ^ttt— + • • • = (D 



du a y ' dD iU a K ^ duy 

by solving the linearized symmetry condition, 

X(E/3(L))=0 on solutions of (^2), P = l,...,q. 

The solution is a g-tuplc Q(x, [u]) = (Q , . . . , Q q ), which is called the characteristic 
of the symmetries generated by X. The set of characteristics of all generalized 
symmetries of the Euler-Lagrange equations (|2.2j) is a Lie algebra, whose Lie bracket 
is 

[Q 1) Q 2 ] = X 1 (Q 2 )-X 2 (Q 1 ). (2.3) 
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Variational symmetries (in the broadest sense) are generalized symmetries whose 
infinitesimal generators satisfy the additional condition 

X(L) = A(^(x,[u])), (2.4) 

for some functions Pq , . . . , Pq. This condition ensures that the variational problem 
is invariant under such symmetries. The set of characteristics of all variational 
symmetries is a Lie algebra with the Lie bracket (|2.3p . Integrating (|2.4[) by parts 
yields an expression of the form 

Q a E a (L) = A(P 4 (x,[u])). (2.5) 

Consequently, each variational symmetry characteristic Q yields a conservation law 

A:(-P l (x, [u])) = on solutions of (123). (2.6) 

This is the local version of Noether's first (and best-known) theorem on variational 
symmetries. 

3. A local approach to Noether's Second Theorem 

Olver (1993) states and proves Noether's Second Theorem in the following form 
(using our notation). 

Theorem 3.1. The variational problem \2. 1\) admits an infinite- dimensional group 
of variational symmetries whose characteristics Q(x, [u;g]) depend on an arbitrary 
function g(x) (and its derivatives) if and only if there exist differential operators 
V 1 , . . . , V q , not all zero, such that 

D 1 E 1 (L) + --- + 2> <Z E 9 (I,) =0 (3.1) 

for all x, u. 

The theorem as stated is restricted to variational problems whose variational 
symmetries do form a group; such problems are uncommon. However, if 'group of 
variational symmetries' is replaced by 'Lie algebra of variational symmetry gener- 
ators,' the theorem applies to all variational problems (|2.ip . We shall use the more 
general wording henceforth. 

The 'if part of Olver's proof is straightforward. Multiply (|3.1j) by an arbitrary 
function g(x) and integrate by parts to get an expression of the form (|2.5I) . where 

Q«(x, [u;g]) = (P«) t (g). 

Then reverse the steps that lead from (|2.4j) to (|2.5[) to show that Q is the char- 
acteristic of a variational symmetry for each g(x). Consequently, if there are no 
characteristics of variational symmetries that depend on a completely arbitrary 
function of x, no relation of the form (|3.ip exists. 

The 'only if part of the proof is quite long, so we do not reproduce it here. 
Instead, we use the following constructive proof, which is much simpler. The key is 
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to regard g as an additional dependent variable and to apply the Euler-Lagrange 
operator E g to (|2.5|) . giving 

E g {Q a (x,[u;g])E Q (L)} =0. (3.2) 

This is the required differential relation between the Euler-Lagrange equations. 
Explicitly, it amounts to 

(-D) J ^ gtt( ; g | U;gl) E a (L)),0. (3.3) 

If Q is linear in g and its derivatives, this relation is independent of g. Otherwise, 
one can linearize Q by taking its Frechet derivative with respect to g, as explained 
by Olver; we discuss an alternative approach later in this section. Typically the 
characteristic is independent of high-order derivatives of g and so (|3.3[) is easy to 
determine; the following example illustrates this. 

[Note: in each example we use whatever notation is standard, so that readers 
can see the results in a familiar form, without having to translate from the notation 
that we use for the general theory.] 

Example 1 The interaction of a scalar particle of mass m and charge e with an elec- 
tromagnetic field has a variational formulation with well-known gauge symmetries. The 
independent variables are the standard flat space-time coordinates x = (x° , x 1 , x 2 , x 3 ), 
where x° denotes time. The dependent variables are the complex-valued scalar ip (which 
is the wavefunction for the particle), its complex conjugate, ip* , and the real- valued elec- 
tromagnetic four-potential J 4 M , /j, — 0,1,2,3; indices are raised and lowered using the 
metric n = diag{— 1, 1, 1, 1}. The Lagrangian is 

L = ijF^F"" + (V M </>) ( V M V0* + m 2 ^* (3.4) 

where 

F^iu = A-u,^ A^^, = Dfj, -\- leAp. 

Therefore the Euler-Lagrange equations are 

= Ey,(L) ee - (V M V»* + ™V, 

= E^. (L) = -V M V M V + m 2 i>, 

= E CT (L) ee ieip(V a ip)* -ieip*V rT ip + na a Ff, 

where E CT is obtained by varying A" . The variational symmetries include the gauge sym- 
metries 

ip h-s> exp(— ie\)ip, ip* \-¥ exp(ieA)-(/>*, A" h-> A* + rf a \, a , 

where A is an arbitrary real- valued function of x; the Lagrangian L is invariant under these 
transformations, so XL — 0. The characteristics of the gauge symmetries have components 

= -ietpg, Q* = ie^*g, Q a = rf a g t0l , 

where g is an arbitrary real-valued function of x. Consequently 

EglQ^ (L) + Q<"* (L) + Q CT E CT (£,)} = -ie^E* (L) + ie^'E*. (L) - D a (^"E, (L)) , 
and so Noether's Second Theorem yields the differential relation 
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Our proof of the 'only if part of Noether's Second Theorem is essentially a 
local version of Noether's original proof (see Noether, 1918). Like Noether's proof, 
it extends immediately to variational symmetries whose characteristics depend on 
R independent arbitrary functions g = (g 1 (x), . . . , g R (x)) and their derivatives. 
This gives R differential relations between the Euler-Lagrange equations: 

E r {g a (x,[u ;g ])E Q (i)} = (-D) J ^Q a (^>;g]) Ea(L) ^| =0i r = i,... )jR . 

(3.5) 

A useful observation is that this set of relations is invariant under any locally 
invertible change of the arbitrary functions. For suppose that each g r is a locally 
invertible function of a new set of arbitrary functions h p , p = 1, . . . , R. Then the 
identity 

E hP {Q«(x,[u;g(h)])E Q (L)} = (E r {Q<*(x, [u; g])E Q (Z)}) (3.6) 

shows that the set of relations 

E h p{Q Q (x,[u;g(h)])E Q (L)} = 0, p=l,...,R, 

is obtained by taking an invertible linear combination of the original relations (|3.5[) . 
So if there exists a locally invertible change of variables such that Q is linear in h, 
(|3.5I) is equivalent to a set of relations that is independent of any arbitrary functions. 

It is convenient to regard Q3.5P as 'new' Euler-Lagrange equations that arise 
when g is varied in the functional 

£[u;g] = j L(x,[u;g])dx, (3.7) 

where 

L(x, [u; g]) - Q«(x, [u; g])E (L(x, [u])). (3.8) 

To prove that the relations (13.51) are independent, suppose that the converse is true; 
then there exist differential operators V r such that 

^Egi (£) + ... + V R V gR (L) = 0. 

By the 'if part of Noether's Second Theorem, there exist variational symmetries 
of (|3.7p whose generator is of the form 

v 1 " 



'{(*T(fi}4- 



where g is an arbitrary function. In other words, there exist functions P'(x, [u; g; g]) 
such that 

1{q q (x, [u;g])}E a (L(x, [u])) - Jt(L) = D^P 1 ). 
Consequently, for arbitrary g and g, 



Q Q (x, [u;g;g])=x{Q Q (x, [u;g])} 
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is a characteristic of variational symmetries for the original variational problem 
(12. lj) . This implies that the set of variational symmetries for (|2.1j) depends on 
more than R independent arbitrary functions, contradicting our original assump- 
tion. Hence the relations (13.51) cannot be dependent. 

Conversely, given R independent relations between the Euler-Lagrange equa- 
tions, Olver's 'if proof applied to each one shows that there are characteristics of 
variational symmetries that depend on R independent arbitrary functions, namely 

Q Q (x,[u;g]) = (P«)V). (3.9) 
Thus we arrive at the general form of Noether's Second Theorem. 

Theorem 3.2. The variational problem \2.1)) admits an infinite- dimensional Lie 
algebra oj variational symmetry generators whose characteristics Q(x, [u; g]) depend 
on R independent arbitrary functions g = (g 1 (x), . . . , g fi (x)) and their derivatives if 
and only if there exist differential operators T>" that yield R independent differential 
relations, 

V?E a (L) = 0, r=l,...,R, (3.10) 

between the Euler-Lagrange eguations. Given the relations h3.10\) , the corresponding 
characteristics are h3. 9\) Conversely, given the characteristics, the corresponding 
differential relations are h3.5\) . 



4. Constrained variational symmetries 

Constraints on arbitrary functions in Q arise from the linearized symmetry condi- 
tion for the Euler-Lagrange equation, coupled with the (linear) requirement that 
the symmetries are variational. Therefore we now suppose that the characteristics 
of variational symmetries of (|2.ip depend on functions g that are subject to S linear 
differential constraints, 

3> sr (g r ) = 0, s = l,...,S, (4.1) 

where each 3> sr is a differential operator. We assume that this set of constraints 
is complete, in the sense that it yields no additional integrability conditions. The 
constraints can be incorporated into the Lagrangian (|3.8[) . yielding 

L(x, [u;g]) = Q«(x, [u;g])E Q (L(x, [u])) - ^ s ,(g r ), (4-2) 

where v , . . . , v s are Lagrange multipliers. Taking variations of (|3.7I) with respect 
to g, we obtain 

(-D)j | ^Q a (x>;g]) e q (L)^ ee (0 flr )V), r = l,...,R. (4.3) 

If S < R, one may be able to eliminate the undetermined functions v s from 
(|4.3p . as will be illustrated by the next example. The elimination of a particular 
set S of dependent variables from a differential system may be accomplished us- 
ing algorithms in differential algebra. Hubert (2000) and Hubert (2003) provide a 
careful and extensive summary of the state of the art to date. Differential elimi- 
nation algorithms have been implemented in the Maple library packages dif f alg 
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and r if simp. These algorithms come with a 'certificate'; provided that the correct 
elimination term ordering is used, the elimination ideal is generated by those equa- 
tions £ in the output that are independent of the variables in S. This means that 
any equation that is algebraically or differentially derived from the input set, and is 
free of the elements of S, can be derived from £ . If £ is empty then the elimination 
ideal itself is empty; in this case, none of the equations that can be derived from 
the input set are independent of all variables in S. 

Incidentally, the elimination algorithm can be used to test directly whether or 
not an Euler-Lagrange system has a differential relation such as (|3.10[) : simply 
replace {E a {L) = 0} with {E a {L) - f a = 0} and eliminate S = {«"}. The result- 
ing equations for {/"} codify the ideal of differential relations on the expressions 
{E a (L)}. 

There are subtleties for nonlinear systems; indeed, the proof that the output 
of the algorithms has the stated properties is one of the most subtle and complex 
problems of differential algebra. 

Example 2 Suppose that Q depends on two functions g 1 (x 1 ,x 2 ) and g 2 (x 1 ,x 2 ) that 
are subject to the single constraint 

g\ + g 2 2 = (4.4) 

(This type of constraint occurs when the variational symmetries include all area-preserving 
diffcomorphisms of the plane.) Then (|4.3|l amounts to 

,_ D)j (^™ E „ (i) ) . _„,„ 

By eliminating the Lagrange multiplier u, we obtain a single differential relation between 
the Euler-Lagrange equations: 

{ . Lh)[ _ Dlj ( Wx,[u;g]) \ +/?i( _ D)j / ^(x> ;g ]) EQ(L) ^ (4 _ 5) 



This is precisely the same constraint as would be obtained by using the local solution of 

El, 

1 2 
g = g,2> g = ~g,l 

(where g is entirely arbitrary) in Noether's Second Theorem. However the Lagrange mul- 
tiplier method has the advantage of working even when (|4.1[l cannot be solved explicitly. 

When the system of constraints (14.1[) does not allow the multipliers v s to be 
eliminated, the most that one can do is to construct conservation laws of the Euler- 
Lagrange equations, as follows. First determine a particular solution of (|4.3I) for u s ; 
this can be done by inspection (which is fast), or by the homotopy method (which 
is systematic). The resulting solution is then substituted into 

v s % r {g r ) ~ g r (^W)V) = when Ei(L) = -.- = E a (L)=0, (4.6) 

yielding a set of conservation laws that depends on g. This set is independent of 
which particular solution is used. It can be shown that if Q is homogeneous with 
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respect to g, the conservation law (|4.6[) is a multiple of the one that is obtained 
from Noether's (first) Theorem. 

Example 3 As we discussed in the Introduction, the Lagrangian 

L=\{{u, x f-{u, t f) 

for the one-dimensional scalar wave equation, 

u,tt — u.xx = 0, (4.7) 

admits variational symmetries with partly-constrained characteristics. Indeed, the most 
general variational symmetry characteristic is Q — g (x,t) + g 2 (x,t), where 



So the modified Lagrangian (|4.2|l is 

L = {u, tt -w^Xg 1 +g 2 ) - ^(g^ + g^i) - ^(g^ -g 2 t )- 
Taking variations with respect to g = (g^g 2 ) gives the identities 

U.tt ~ U <xx = ~{y] x + v] t ), U <tt ~ U.xx = - ^t)' 

which are satisfied by 

v 1 = u iX —u y t, v 2 — u tX + u tt - 
Then the conservation laws (j4.6[) amount to 

{(g 1 +g 2 )«, 3! -(g 1 -g 2 )«,t} fX + {(g 1 -g 2 )«^-(g 1 + g a K*} )t =0 when g2D holds. 
There are no trivial conservation laws with Q nonzero. 

Example 4 For a more substantial example that produces interesting conservation laws, 
consider the two-dimensional shallow water equations from Lagrangian fluid mechanics. 
Each fluid particle is labelled by its position (a 1 , a 2 ) at some reference time, t = 0, and 
(x,y) = {x(a} , a 2 , t), y(a} , a 2 , t)) is the position of the fluid particle at time t. We set 
x = x t t and y = y,t', from here on only the derivatives with respect to the label variables 
a 1 are denoted by subscripts. Salmon (1983) introduced a Lagrangian for the shallow water 
equations; with constant Coriolis parameter /, it amounts to 

y]) = g ( ±2 + y 2 ) + fay ~ \s h > 

where 

h= , 

X,lV,2 — X,2V,1 

is the fluid depth and g is the constant of gravity. To translate between the Eulerian and 
Lagrangian viewpoints, we need the identities 

8 , / 8 8 \ 8 , f 8 8 



dx \' 2 da 1 da 2 J ' dy \ '"da 1 ' ' L da 2 

Calculating the Euler-Lagrange system yields the shallow water equations 

Bh Bh 
E x (L) = -x + fy~g—^0, E y {L) = -y - fx - g — = 0. (4.8) 
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These are supplemented by the continuity equation, 

which is not an Euler-Lagrange equation in Salmon's formulation. The variational sym- 
metries include the well-known particle relabelling symmetries; these are arbitrary area- 
preserving diffeomorphisms of label space, whose characteristics are (Bila, Mansfield & 
Clarkson, 2006) 

Q X = x ,ig\ Q V = V,ig\ where g 1 = g 2 = g^ + g 2 2 = 0. 

Taking the three constraints on the functions g = (g^g 2 ) into account, the modified 
Lagrangian is 

L = x^gE x {L) + y^gEy(L) - ^g - v 2 g - t/ 3 (g\ + g 2 2 ) . 
Taking variations with respect to g, we obtain 

xjE x (L) + y,iE v {L) = -v - v% i = 1,2, 

which is satisfied by 

v 1 = xx,i + yy,i +fxy,i, 
v 2 = xx, 2 + yy,2 + fxy,2, 
v A = -|(i 2 +y' 2 ) - fxy + gh. 

With these solutions i/, the resulting conservation laws (14. 6p arc 

A(gV +gV) + (gV 3 ),i + (gV), 2 = when gU) holds. 

In particular, taking (g^g 2 ) to be (1,0) and (0, 1) in turn yields 

u 1 + u\ = 0, i> 2 + v\ = when holds, 



which were discovered in a (somewhat complicated) multisymplectic setting by Hydon 
(2005). Differentiating the first of these with respect to a 2 , the second with respect to a 1 
and eliminating ^ 3 12 , we obtain the well-known conservation law for potential vorticity, 

D * \ T (lT ~ IT + f) \ = when <E3 holds. 



h \dx dy 

So in this example, our extension of Noether's Second Theorem has led quickly and easily 
to some fundamental conservation laws. Conservation of potential vorticity also arises 
directly from Noether's Theorem, but considerable effort is needed to obtain it (see Bila, 
Mansfield & Clarkson, 2006). 



5. Prom continuous to discrete: Noether's Second Theorem 

For difference equations, the dependent variables are again u = (u , . . . , it 9 ), but 
now the independent variables are n = (n , . . . ,n p ), where each n % is an integer. 
The forward shift operators, Si, are defined by 

Si-.*? -mi + Si, $/(n) = /(^n), 
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for all functions /; here <K is the Kronecker delta. The adjoint of S, is the backward 
shift (S'. i ) t = S7 1 : n j h4 n j - 5{. We also use the multiple shifts, Sj = Sf • • • S J P P , 
and the forward difference operators, 

A = Si - id, 

where id denotes the identity map. (Here and in what follows, we use notation 
that is similar to the continuous case in order to highlight the similarities. To avoid 
confusion, we use a tilde to distinguish an operator from its continuous counterpart.) 
For difference variational problems, the action is of the form 



C[u] = $>(n, [u]), 



where now [u] denotes u(n) and finitely many of its shifts. Then the Euler-Lagrange 
equations are 

E Q (i) = S_j |^-|=0, a = l,...,q. (5.1) 

Generalized symmetries of the Euler-Lagrange equations have generators of the 
form 

X = S J (Q«(n I [u]))^ 

that satisfy the linearized symmetry condition, 

x(E a {L)) = when fl5T]) holds. 

Here Q Q (n, [u]), a — 1, ... ,q, are the components of the characteristic, Q(n, [u]). 
The symmetries are variational if they do not change the action, that is, if there 
exist functions Po(n, [u]) such that 

X(L)=DiI$(n,[u]), (5.2) 

Summing (|5.2[) by parts yields the following condition, which is analogous to (12. 5p : 

Q a E a (L)4,P'(n,[u]), (5.3) 

for some functions P\ So Noether's Theorem applies just as in the continuous case: 
every variational symmetry characteristic yields a conservation law of the Euler- 
Lagrange equations, namely 

D,P(n, [u]) = 0, when dHH) holds. 

With the above notation, the analogue of Noether's Second Theorem for the 
difference calculus of variations is as follows. 

Theorem 5.1. The difference variational problem <5 u £[u] = admits an infinite- 
dimensional Lie algebra of variational symmetry generators whose characteristics 
Q(n, [u;g]) depend on R independent arbitrary functions g = (g 1 (n), . . . , g fl (n)) 
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and their derivatives if and only if there exist difference operators T>" that yield R 
independent difference relations, 

V?E a (L)=0, r = l,...,R, (5.4) 

between the Euler-Lagrange equations. Given the relations |5.^| ), the corresponding 
characteristics are 

Q Q (n,[u;g]) = (P«)V). (5.5) 
Conversely, given the characteristics, the corresponding difference relations are 

E g .{Q«(x, [u;g])E Q (Z)} = S_j ( ^ Q g ( si" ;gD ^ (L) ) = °' (5 ' 6) 
and consequently 

m=<S VWx,[u; g] ) 



The proof completely replicates our proof of the continuous version, so we shall 
omit it. The identity (|3.6p holds for difference equations, so the set of relations is 
invariant under locally invertible changes of the arbitrary functions. 

Example 5 For the interaction of a scalar particle with an electromagnetic field (see 
Example 1), Christiansen & Halvorsen (2011) discovered a finite difference approxima- 
tion that preserves the gauge symmetries as variational symmetries. The mesh x(n) is 
uniformly-spaced in each direction, with step lengths 

ft" = (S„- id) (x M (n)), = 0,..., 3. 

It is useful to introduce the scaled forward difference operators 

D„ = D li /h fl , = 0,..., 3, 



the adjoint of the (scaled) forward difference operator is the negative of the (scaled) 
backward difference operator. Let ?/>(n) denote the approximation to the wavefunction tp 
at x(n). Christiansen and Halvorsen used a Yee discretization to approximate on the 
edge that connects the points x(n) and S M x(n); we denote this approximation by A M (n). 
As before, indices are raised and lowered using the metric n = diag{ — 1, 1, 1, 1}. Up to a 
sign, the Lagrangian for the scheme is 

L = jF^F^ + (v M ^) (v"#) * + m 2 ^* (5.7) 
where, for each n, v = 0, . . . , 3, 



V M = tJj {S M - exp(- ieh^A^id} = + ±- {l - exp(- ie^l M )} 



id. 
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Therefore the Euler-Lagrange equations that are obtained by varying ip, ip* and A" are 
= E ? .(L) = (V^V^ + m 2 ^, 

= E(j (Z/) ee ieexp(- ie/i <T l (T )^(V CT ^)* - ie exp(ie/i CT 2 CT )^* V CT -0 - ij^c'r". 
Note that 

V* = 5* + ^ { 1 - exp ( - ieft" A^) } id 

tends to — V£ as W approaches zero, so the scheme has the correct limiting behaviour. 

The above discretization preserves the gauge symmetries as variational symmetries. 
Explicitly, these are 

ip h-s> exp(— ieA)^>, ip* h-» exp(ieA)f/>*, A a h-s- A" + rf a D a \, 

where now A is an arbitrary real-valued function of n. As in the continuous case, the 
Lagrangian is invariant, so XL = 0. The characteristics of the gauge symmetries have 
components 

= -ieipg, = ie^*g, Cf = rf a D a g, 

where g is an arbitrary real-valued function of n. Consequently the discrete version of 
Noether's Second Theorem immediately yields the difference relation 

-ie^E^(£) + ie^*E f .(L)+S* (», aa E„(L)) ee 0, 

which Christensen and Halvorsen obtained with the aid of Noether's first theorem. 

6. Variational symmetries with difference constraints 

If the functions g are subject to a (complete) set of S linear difference constraints, 

@srg r =0, S = 1,...,S, (6.1) 

we incorporate these using the same approach as for differential equations. First 
form the new Lagrangian, 

L(n, [u;g]) = Q a (n, [u; g])E Q (L(n, [u])) - v s § sr (g r ), (6.2) 

where v , . . . , V s are Lagrange multipliers. Then take variations with respect to g 
to obtain 

S-a ( dQC ^ g]) %.{L)) = feV), r = l R (6.3) 

Just as in the continuous case, these relations yield conservation laws of the Euler- 
Lagrange equations (15.11) . namely 

u s 3> sr g r - g r (3>sr)H" s ) = on solutions of (6.4) 

It may also be possible to eliminate the Lagrange multipliers from (|6.3p to obtain 
difference relations between the Euler-Lagrange equations. 
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Example 6 The lattice Korteweg-de Vries (KdV) equation is 

U 1A -U , -c( _L-_iL) =0, (6.5) 



t/i,o U , 



where Uij denotes the value of the dependent variable U at the point (n 1 + i, n 2 + j); 
see Grammaticos ei ai. (1991) for details. This is not an Euler-Lagrange equation as it 
stands, but it can be turned into one by introducing a potential, u, such that 

Then ()6.5|l amounts to 

E U (L) = w ll0 - uo,i + w-i.o - «o,-i - c 

where 



Ui.-i — uo,o uq,o — U-1,1 



L[u] = 1*0,0(1*1,0 - i*o,i) + cln(wi,o - 1*0,1). 
Clearly, the variable U is unaffected by gauge transformations 

Uij t-¥ iti,j + Ai,j, where Ai,o = Ao,i, 

so these are symmetries of the Euler-Lagrange equation. Their characteristic is Q = g , 
where g x — g x = 0; therefore 

X ( L ) = go,o( w M - «o,i) = L>i(g_ lj0 i*o,o) + g ,-i no,o), 

and so the gauge symmetries are variational. Now vary g for the new Lagrangian 
L[u; g] = QE u (L[u]) - iA>,o(g 1>0 - g ,i) 

to obtain the relation 

E u (Lki]) = i/_i,o - fo,-i- 

One obvious solution is 

fo.o = ito.o - ?ti.i H ; (6.6) 

1*1,0 — 1*0,1 

the corresponding conservation laws (|6.4|1 are 



Di(g V-i,o) + D 2 (- g fo,-i) = when E U (L) = 0. 



Note that by setting uo,o = in (|6.6|l . we obtain the well-known discrete potential KdV 
equation, 

(t*i,i — wo,o)(«i,o — wo,i) = c. (6.7) 

So the constraint on the gauge symmetry provides a link between the lattice KdV equa- 
tion and its related potential form. This observation applies similarly to many integrable 
difference equations. 



7. Conclusion 

We have shown there is a complete correspondence between Noether's Second The- 
orem for differential and difference equations, and that the same is true where 
characteristics of variational symmetries depend on partly-constrained functions. 
Our approach yields, mutatis mutandis, the corresponding results for differential- 
difference equations. 
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